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The reaction mechanism of projectile fragmentation at intermediate energies has been investi- 
gated observing the target dependence of the production cross sections of very neutron-rich nuclei. 
Measurement of longitudinal momentum distributions of projectile-like fragments within a wide 
range of fragment mass and its charge was performed using a hundred- Me V/n *°Ar beam incident 
on Be and Ta targets. By measurement of fragment momentum distribution, a parabolic mass de- 
pendence of momentum peak shift was observed in the results of both targets, and a phenomenon 
of light-fragment acceleration was found only in the Be-target data. The analysis of production 
cross sections revealed an obvious enhancement of the target dependence except target size effect 
when the neutron excess is increased. This result implies the breakdown of factorization (BOF) of 
production cross sections for very neutron-rich nuclei near the drip line. 

PACS numbers: 25.70.Mn, 25.70.-z, 27.20.-fn, 27.30.-|-t, 27.40.4-z 



I. INTRODUCTION 

During these two decades, unstable nuclear physics 
has become one of the most interesting fields of nuclear 
physics. Nuclear fragmentation of heavy-ion beams is 
utilized for producing secondary beams of unstable nuclei 
far from /?-stabihty. For designing experiments with the 
secondary beams, a good knowledge of production cross 
sections is essential. To deduce the production cross sec- 
tions relevantly, an empirical parametrization of fragmen- 
tation cross sections (EFAX) is widely used in simula- 
tion programs for projectile- fragment separators [H, HjQ. 
For consideration of reaction mechanism, fragmentation 
models based on abrasion-ablation scheme (AA models) 
are often applied to estimate the production cross sec- 
tions. 

With the recent development of heavy-ion accelerators, 
the number of accessible nuclei lying far from /3-stability 
line has been increasing. An example has been shown by 
the experimental findings of particle stability of ''"'^F [j], 
3i'34Ne ^, ^6(], 37Na ^6,], 37,38Mg I, g, 40,41^1 and 
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Si [6[, as the most neutron-rich nuclei so far identified 
for the Z=9— 14 elements. For the new isotope-search 
experiment, the accuracy to predict the production cross 
section of neutron-rich nuclei near the drip line is very 
important to discuss the particle stability of them. How- 
ever, the predictive power of the EFAX parametrization 
and AA models is not strong enough for specific very 
neutron-rich nuclei 8]. For instance, although tantalum 
is often experimentially used as a production target in 
order to earn the better yield of these nuclei, the target 
dependence of production cross sections is not taken into 
account in these models. Instead, we first determined the 
production cross sections for the observed isotopes, which 
were then used to estimate the production cross sections 
and the expected yields for the non-observed isotopes. 
Therefore, to enhance predictive powers for production 
cross sections, deeply understanding reaction mechanism 
as well as observing systematic behaviors of production 
cross sections are necessary. 

The validity of the EFAX parametrization and AA 
models has been mainly verified for medium and heavy 
mass fragments via multi-GeV high-energy fragmenta- 
tion reactions. The EFAX formula has been also used 
for the intermediate-energy experiments (several tens A 
MeV) since the formula can reproduce reasonably well 
even at intermediate enegies the production cross sec- 
tions of stable and unstable nuclei near the /3-stability 
line. A target dependence of fragment-production cross 
sections in relativistic heavy-ion collisions is limited only 
by taking account of the nuclear-size effect [9,]. However, 
the recent experiments show that the production yields 
of nuclei far from stability line are quite different from 
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the prediction of EPAX formula, and strongly dependent 
on the N I Z ratio of the target [lOj. It is of great interest 
whether the cross sections measured with different tar- 
gets factorize in projectile fragmentation at intermediate 
energies. 

In order to confirm the validity of factorization, a care- 
ful measurement of fragment momentum distribution is 
necessary for precise determination of the production 
cross sections. At relativistic energies, the shape of mo- 
mentum distribution for an isotope was found to be a 
Gaussian function 11] and the width was well understood 
with a statistical model [l2|. On the other hand, the 
momentum distribution of fragments produced at inter- 
mediate energies has an asymmetric shape with a tail at 
low momentum side. A theoretical attempt was made to 
reproduce the asymmetric shape, by taking into account 
nucleon flows between projectile and target in the colli- 
sion time |T^. In this model, stochastic nucleon trans- 
fers using a Monte Carlo method and sequential evapora- 
tion were taken into account. This model can reproduce 
the low-momentuni-side tails, however, due to a large 
friction force, the whole of predicted distributions are 
shifted toward the low momentum side much larger than 
the observed. This discrepancy of momentum distribu- 
tions affords a large ambiguity to evaluate the production 
cross section from the measured yield of fragments. In 
addition, the measurement for very neutron-rich nuclei 
has been performed using very thick targets to earn the 
yields, so far. As distortions of momentum distributions 
due to target thickness, the momentum distributions of 
fragments cannot be obtained from these data clearly. 
Therefore, the measurement of momentum distribution 
of very neutron-rich nuclei has become important. 

In the present work, we focus on the target depen- 
dence of momentum distribution of projectile-like frag- 
ment (PLF) produced by nuclear fragmentation reactions 
at an intermediate energy. To investigate the target de- 
pendence of the production cross sections systematically, 
we used two production targets of Be and Ta. To avoid 
distortions of momentum distributions due to the target 
thickness, we prepared relatively thin targets. We per- 
formed the experiment with the RIKEN-RIPS to elimi- 
nate the primary beam and to collect the projectile-like 
fragments. The data of this experiment were taken in a 
wide range for fragment mass including very neutron-rich 
side {N / Z J « 3) toward the neutron drip-line and light 
mass {Af > 3), with a good statistics for momentum- 
distribution tails. 

In the following, we first describe the experimental 
setup and procedure in Sec. II. In Sec. Ill, the analy- 
sis of the data is described. In Sec. IV, the experimental 
results are described, where the observed momentum dis- 
tributions of fragments, the momentum peak shift, the 
high- and low-momentum side widths, and the produc- 
tion cross sections are presented in terms of target depen- 
dence. Based on the results, we discuss the prefragment 
production mechanism in projectile fragmentation reac- 
tions. 



II. EXPERIMENTAL SETUP AND 
PROCEDURE 

The projectile fragmentation experiment using a 40 Ar 
beam was performed at the RIKEN Accelerator Research 
Facility. The measurement of momentum distributions 
of projectile- like fragments was carried out with the pro- 
jectile fragment separator, RIPS Q^. The 40^j.i7+ ]-,gaj^ 
accelerated by the ring cyclotron with energies up to 90 A 
MeV and 94A MeV irradiated a 95-mg/cm^ thick ^Be 
target and a 17-mg/cm^ thick '"^*Ta target, respectively. 

The primary beam intensity was monitored for nor- 
malization of fragment yields to obtain the momentum 
distribution. A plastic telescope consisted of three plastic 
scintillators with a size of 50x50x0.5 mm'^ was placed at 
a backward angle of 135 degrees and at a distance of 0.5 
m from the target position. The plastic detectors counted 
yield rates of light particles produced with nuclear reac- 
tions at the production target. The three photomultiplier 
tubes (PMTs) were mounted on the plastic scintillators 
one-by-one. The counting rate of triple coincidence was 
used to monitor the primary beam. The primary beam 
intensity was measured with the indirect method of plas- 
tic counters calibrated by a Faraday cup. By changing 
the primary-beam intensity which ranges from 

10^4 to 

1.0 of the full beam intensity, the calibration data were 
taken. The systematic error of the beam monitor was 
estimated to be 7%. 

The RIPS was used as a doubly achromatic spectrom- 
eter. Projectile fragments produced at the production 
target and emitted at 0° were collected and transported 
to a double achromatic focal plane (F2). The momentum 
acceptance was set to be Ap/p=l% at a momentum dis- 
persive focal plane (Fl) where left and right slits formed 
the rectangle window of momentum acceptance. The an- 
gular acceptance was set with a square window formed 
by four slits (upper, lower, left and right), which were 
placed at the behind of the production target. The 9 
and (j) angular acceptances were 25 mrad, which is nar- 
rower than the width of angular distribution of fragments 
in r.m.s. We use the constant value of 6.25x10^^ msr as 
the solid angle of Ail. 

Momentum distributions of fragments were measured 
at 23 settings of magnetic rigidity {Bp) over a range of 
2.520-4.068 Tm using ^Be target. In the case of ^^^Ta 
target, the measurement of momentum distribution for 
each fragment was performed at 31 magnet settings as 
same as Be case. When the magnetic field was changed 
for each run, the F2 image of secondary beam was mea- 
sured by means of a parallel-plate avalanche counter 
(PPAC) [l^ to confirm transmission, and the x-position 
of the beam at F2 was corrected to center precisely by 
tuning the D2 field. The beam position was monitored 
with an accuracy of 1 mm. The systematic error of mag- 
netic rigidity was about 3x10"^ from the ratio of 1 mm 
to 3.6 m. The difference between Dl and D2 magnetic 
fields was less than 0.05%. 

The identification of fragments was carried out event- 
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by-event by means of measurement of time of flight 
(TOF) and energy deposit (AE) for each fixed Bp run 
with the 1% momentum sht. According to an estimation 
of charge state distribution ^16*1 , aU fragments in flight are 
fuUy striped {Q=Z). Under this assumption, the parti- 
cle identiflcation can be performed on the basis of the 
relations: 



A 

Bp (X -13 



(1) 
(2) 
(3) 



where A and Z are mass number and atomic number, 
respectively. 

The detectors of two 0.5-mm thick surface-barrier-type 
silicon counters (SSDl, SSD2) and a 0.5-mm thick plastic 
scintillation counter (PL) were installed for the AE and 
TOF measurement at F2. To use two silicon detectors 
allowed us to deduce Z number of the fragment indepen- 
dently, and to take correlations between them to achieve 
good S /N ratio. 

Both SSDl and SSD2 have 48 mmx48 mm sensitive 
area which is wide enough to accept all particles reaching 
to F2 where FWHM of the beam profile is 6 mm. Two 
photomultiplicr tubes (PMTs) were mounted on both 
sides of PL (Left and Right). A timing of PL was de- 
termined with an average of the signals from the Left 
and Right. 

The TOF of each fragment over 21.3-m fiight path be- 
tween the production target and the F2 was determined 
from the difference of timing signals between RF signal 
of the cyclotron and the PL timing. The TOF resolution 
was measured with a faint beam to be 0.27 ns (r.m.s.), 
which included the timing jitter of RF signal ('^0.09 ns). 
Thus, the intrinsic resolution for PMT was estimated to 
be 0.18 ns. 



III. DATA ANALYSIS 

In this section we describe the procedures of the data 
analysis, fitting of momentum distributions, and evalua- 
tion of the production cross sections. 



A. Particle identification 

Figure [1] shows a two-dimensional plot versus AE in 
SSDl for one Bp setting using the Be target. A rejection 
of the background events was carried using a correlation 
gate between SSDl and SSD2. By use of the 3-cr gate 
by two SSDs, we achieved the particle identification of 
fragments with low background events. Fragment yields 
were obtained by counting the isotopes from the particle 
identiflcation. Figure [5] shows the Z- and ^/Z-projection 
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FIG. 1: Particle identification in the dE-TOF plane at Bp ■ 
3.708 Tm with the Be target. 



of the particle identiflcation at _B/9=2.523 Tm setting us- 
ing the Be target. The counting gate of isotopes was a 
rectangle region with ±3-cr of the resolution (^(r.m.s.) for 
Z and A/Z. 

In the data acquired using the Be target, the analyzed 
isotopes were ^-SLi, ^-i^g^^ lo-isg^ 11-18^, 13-21^^ 

15-240, 17-27p^ 19-29Ne, 21-32Na, ^^-^^Mg, 25-36 Al, 
27-38si, 29-39p^ 33-38s, 36-39^1^ ^nd ^^kr. Most of 

these isotopes are neutron-rich ones. It should be 
noted that the 3'5Al,37,38gj^ 38,39p jgotopes ana- 

lyzed have neutron numbers larger than the projectile 
(N > 23), which are produced through neutron pick- 
up process. We analyzed from the Ta-target data for 
6-8Li, 9-iiBe, lo-i^B, ii-i^C, 13-19N, i5-2iO, 17-24F, 

19-27Ne, 21-29Na, 23-3lMg, 24-34^1, 26-34gi^ 29-36p^ 

30-37S, 33-37ci^ 35-39 Ar, and 37-40^. The potassium 
isotopes should be produced by the reaction with pro- 
ton pick-up process. 

In order to obtain the doubly-differential cross section 
from each fragment yield, we estimated the transmission 
between FO to F2, and the reaction loss in the detec- 
tors. A Monte Carlo simulation by MOCADI [2] was 
performed under the realistic condition of RIPS using a 
reference beam of '^'^Ar. The value of transmission was 
obtained to be 95.3±0.3%. The nuclear reaction loss of 
fragment in the detectors was evaluated with reaction 
cross sections calculated by a simple geometrical model. 
The reaction loss in the detectors was estimated to be 
less than 0.8%. After all, the systematic error was ±9% 
for evaluation of the fragment cross sections. 



4 



10 



4 _ 



10 



B 

o 

loi 



3 _ 



10 



I 



Ms 




6 8 10 12 14 

Charge, Z 



10- 



X2 



O 



10 



2 _ 



10^ 




U 10 1 -^*A1 



FIG. 2: (a) ^-projection spectrum for A/ Z—2it0.3 and 
(b) A/Z-projection spectrum for Z=13.0±0.5 at Bp— 2.523 
Tm (Be target). The arrow in (a) indicates the lack of ^Be, 
which is known to be particle unbound. 



B. Fitting Procedure 

Momentum distributions of fragments have informa- 
tion for understanding reaction mechanisms. At rela- 
tivistic energies, the projectile fragments have symmet- 
ric momentum distributions fitted to a Gaussian form, of 
which width has been discussed with respect to the Fermi 
motion of nucleons or temperatures of pre-fragments . 

The results obtained in this work are different from at 
relativistic energies. Figure [3] shows a typical momen- 
tum distribution of this experiment. By comparison of 
the momentum distribution fitted with a Gaussian func- 
tion (dotted curve), asymmetric feature of the distribu- 
tion is clearly observed. The momentum distribution of 
projectile-like fragments produced at intermediate ener- 
gies are generally asymmetric with a tail on the low mo- 
mentum side ^17.]. 

To deduce the most probable momentum and width 
from such skewed shapes, the momentum distributions 
have been fitted so far with several kinds of trial func- 
tions [13, m, [3 ■ Since physical models have not been 
established for the low momentum tail, what kind of func- 
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FIG. 3: Typical fragment momentum distribution (dashed 
curve) and fitting results for the momentum distribution of 
^°Mg data. The fitting result with a Gaussian function (dot- 
ted curve) shows clearly asymmetric feature of the experimen- 
tal data. The solid curve indicates a fitting result with the 
asymmetric Gaussian function. 



tions to be used is not unique. To study a systematics of 
low-momentum tail, the following asymmetric function 
with four free parameters is applied for the present data 
fitting. 



1 = 0°) = 



,(4) 



A-e^VK w-2 



where Pq is the most probable peak value of momentum 
in the distribution, gl and an are momentum width in 
low and high energy side, respectively. 

In this fitting procedure, the maximum likelihood 
method was used to treat the small statistics at the tail 
parts of distribution. The result of fitting with the asym- 
metric Gaussian function is shown with the solid curve 
in Fig. El 

By means of this method, we obtained fitting results 
of momentum distribution for all isotopes available in 
our data, though there is another kind of complexity 
in a few case. We found two components in the mo- 
mentum distributions of very light fragments only for 
the Be-target data. Figure [J] shows the momentum dis- 
tributions of the ^°Be and "^^Mg isotopes from the Be- 
target data. The both distributions are scaled as a 
function of velocity (/?). Two arrows in the figure in- 
dicate the velocities of projectile {(3proj ) and center-of- 
mass system (Pcm), respectively. In the distribution of 
30Mg (Fig.[4l[b)), a single component is observed near the 
projectile velocity (Pproj)- On the other hand, the ^"^Be 
distribution (Fig. Hl^a)) shows two components at Pproj 
and Pern- Here, the component around l3proj is defined 
as a high-energy side of peak (HE), and that of f3cm is 
as a low-energy side of peak (LE). We made an attempt 
to fit the data using the asymmetric Gaussian function 
for the HE-component and a Gaussian function for the 
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FIG. 4: Typical fragment momentum distributions for 
(a) ^°Be and (b) ^°Mg produced in Ar+Be reaction. The 
fitting results are also shown with solid curves. Light frag- 
ment of ^"Be has two components (dashes curves) of HE- and 
LE-components, while heavy fragment of '^''Mg has one com- 
ponent corresponding to HE-component. The IfE-component 
for each isotope has a low momentum tail and the symmetric 
parts are shown with dotted curves. 



LE-component. The fitting results are drawn with solid 
curves. With increasing the fragment mass number, the 
LE-component decreases very quickly. At last, no signif- 
icant LE-component has been observed for heavy frag- 
ments like the ^°Mg data. In our data, we have found 
clearly the LE-component for the fragments with A of 
9^12. The LE-component has been observed for light 
fragments in Ar-|-Be reaction, while, the LE-componcnt 
has not been found in the momentum distributions of 
Ar-fTa reaction data in the momentum region where the 
experimental data were taken. When a light fragment 
like ^"B is produced in the Ar-|-Ta reaction, the impact 
parameter is much larger than in the case of the Ar-|-Be 
reaction. Thus, we found the LE-component only for the 
Ar-|-Be system. Since we focus on the projectile fragmen- 
tation reaction, we discuss mainly the HE-componcnt in 
this paper. 



C. Evaluation of cross sections 

The production cross sections of fragments were eval- 
uated with the fitting results of the longitudinal momen- 
tum distributions. The transverse momentum distribu- 
tions are assumed to be a Gaussian function with a width 
of a±, 
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FIG. 5; Overall systematic error for the production cross 
sections (solid line). The dashed line shows error originated 
from ambiguity of a±. 



where an is a parameter of the deflection effect, and 
Apj are mass numbers of a projectile and a fragment, 
respectively. In Ref. the deflection parameters for 
two targets of ^'^Al and ^^'^Au with 92.5 and 117.5A 
MeV ^^O beams. The results from their measurement 
were CT£)=190— 220 MeV/c, found no large target depen- 
dence. Taking account of an energy dependence reported 
in Ref. [Hj, we used ao of 195 MeV/c around 90^ MeV 
for the present experiment. The ambiguity of transverse 
momentum distributions was taken into account as the 
systematic error of the cross sections. The overall sys- 
tematic error for the production cross sections is shown 
in Fig. [3 



IV. RESULTS AND DISCUSSION 

The fitting results with the asymmetric Gaussian func- 
tion are presented in the following. First of all, the fitting 
results of momentum distributions are compared with 
several formulae taken from reaction models. The result 
of high momentum side width can be understood by the 
Goldhaber model. However, we need the further discus- 
sion of the momentum peak shift and low momentum side 
width. Next, we show the result of production cross sec- 
tions. The result reveals the phenomena, breakdown of 
factorization (BOF), for production of very neutron- rich 
nuclei. We discuss the systematics of isotope produc- 
tion cross sections. The charge distribution of the cross 
sections for a fragment mass is obtained from our data, 
and compared with the EFAX formula. Finally, we dis- 
cuss the prefragment production mechanism in projectile 
fragmentation reactions to search for the origin of BOF. 



A. Momentum peak shift 

In the nuclear fragmentation process, a part of the 
kinetic energy of projectile is converted into excitation 
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FIG. 6: Momentum peak shift of fragments produced in (up- 
per) Ar+Be and (lower) Ar+Ta. The primary beam energies 
are drawn with dotted lines, and the kinetic energies at the 
half point of target thickness (projectile energies) are the solid 
lines. See the text for the labels and curves. 



energies of fragments, and the projectile velocity is de- 
creased. This energy loss in projectile fragmentation re- 
action is called as 'momentum peak shift'. The momen- 
tum peak shift is obtained from the difference of the pro- 
jectile velocity and the most probable velocities of frag- 
ments. We present the result of momentum peak shift 
in a unit of energy per nucleon, which is proportional to 
square of velocity. This unit is convenient to discuss the 
kinetic energy consumption in the nuclear reaction. 

Figure [6] shows the momentum peak shift of fragments 
produced in Ar-|-Be and Ar+Ta reactions. The solid lines 
for each isotope are drawn to guide the eye. The primary 
beam energies were corrected with the mean energy loss 
in the production target. The obtained projectile ener- 
gies (solid lines) were taken to be 87.5^ MeV and 93.8A 
MeV, respectively. The values of the most probable en- 
ergies of fragments were also corrected with energy losses 
in the targets. The negative shift of data from the solid 
line of the projectile energy indicates the energy loss by 
nuclear fragmentation process. The systematic error of 
momentum shift was estimated to be 0.9% in the unit of 



energy per nucleon. 

The deviation for a given element chain in the Be tar- 
get is comparable to the systematic error. As seen in 
Fig. [6l the measured momentum shifts for the Be tar- 
get (a), compared with the case of the Ta target (b), 
vary widely. One could understand the deviation as an 
effect of the target thickness in atomic energy-loss pro- 
cess [22!| , because the target thickness of Be is six times 
thicker than that of Ta. 

The bounce of momentum peak shift as a function of 
fragment mass was observed in the results of both targets. 
The momentum peak shift increases when the number of 
removed nucleons AA = Ap — Af is increased up to a 
half of projectile mass {Af >20). On the other hand, 
the momentum shift decreases when the mass loss AA is 
increased beyond 20 {Af <20). In short, we observed the 
maximum of momentum peak shift around Af ^20. 

A phenomenon of fragment acceleration was observed 
in the Be-target data. In Fig.[SJa), the solid line of 87.5^ 
MeV corresponds to the primary beam energy. Veloci- 
ties of the fragments ^'^Li, ^""'^^Be, and ^^.^ larger 
than the projectile velocity. The very light fragments are 
accelerated in the reaction process. We note that the ac- 
celeration phenomena was also observed in collisions of 
238U at 1 A GeV with lead, reported by Enqvist [H]. 
On the other hand, no acceleration phenomenon for all 
fragments Af >8 was observed in the Ta-target data. 

Except the acceleration phenomenon, no significant 
difference between both targets was observed in the 
momentum peak shift for the projectile-like fragments 
Af >2G. 

We found two features in the momentum peak shift. 
First, the maximum of momentum peak shift is observed 
for the fragments around Af —20 in Fig. [6] for both tar- 
gets. Secondly, the acceleration phenomenon was ob- 
served in light fragments only for the Be target. Namely, 
the most probable velocities of light fragments are be- 
yond the projectile velocity. In the following, we discuss 
the features observed in the momentum peak shift. 



1. Parabolic mass dependence of peak shift 

The momentum peak shift has been investigated for 
long time. Many of the reports have shown a linear 
mass dependence of the peak shift for fragments with 
Af > Ap/2. Several formulae are proposed so far to 
reproduce the momentum peak shift for the heavy frag- 
ments produced in peripheral collisions. We compare our 
experimental results for the wide fragment mass range 
with the formulae and make an attempt to introduce a 
new picture to reproduce the parabolic mass dependence. 

In Fig. [51 four formulae are shown with dashed curves 
and the labels as Greiner75, Kaufman82, Borrel85 and 
Morrissey89 [H, [H [H, [H [H, [H, [H . For every for- 
mula, the momentum shift becomes large when the num- 
ber of removed nucleons AA is increased up. This ten- 
dency conflicts the present result in the region oi Af <20. 
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To compensate for the deviation, J. A. Winger et al. [ij ob- 
tained a formula symmetrized mathematically with re- 
spect to Ap/2 (dotted line). However, it still overesti- 
mates the momentum shift around Af ^ Ap/2. After 
all, all of formulae above cannot predict the parabolic 
mass dependence. 

The previous investigations were mainly performed in 
the region of fragments close to the projectile mass. For 
whole fragment mass region observed in this work, a 
parabolic dependence, where the symmetric point is at 
the half mass of projectile, is observed. As Borrel com- 
mented first in Ref. [30|, the symmetric behavior of the 
velocity shift implies another mechanism less costly, com- 
pared with the removal of individual nucleons. The mech- 
anism may be a process that the projectile splits into two 
pieces. Then we make a new formula which calculates the 
momentum shift based on the splitting picture. 

First, we assumed that the projectile splits in two frag- 
ments which are not exited state. By means of the empir- 
ical mass formula of Weizsacker-Bethe M{N,Z), where 
{N,Z) is the neutron and charge number of nuclei [3]| . 
the energy loss in the splitting process is described as 
M{Nf,Zf)+ M{Np - Nf,Zp - Zf)- M{Np,Zp). The 
change of binding energy shows the parabolic mass de- 
pendence of peak shift. However, the energy loss by the 
splitting process gives only several tens percent of the 
value measured as momentum peak shift, quantitatively. 

Next, we deduce a semi-empirical formula in consider- 
ation of the excitation energy of fragments. The splitting 
process consumes the kinetic energy of projectile due to 
the excitation. Here we assume that the energy con- 
sumption is proportional to the number of pairs of nu- 
cleons destroyed in the reaction, where the nucleons are 
acting mutually by long range force. In this picture, the 
number of pairs of nucleons in projectile is counted with 
ApC2=Ap{Ap— l)/2, where the mark of „Cm is the com- 
binatorial which gives the number of ways of choosing 
m out of n. When the projectile splits into the spec- 
tator Af and the participant {Ap — A/), the number of 
pairs of nucleons which decreased by reaction serves as 
Af{Ap — Af ). When one nucleon is removed {Af—Ap—1), 
the number of pairs which decreased is 1, and the 
total energy loss can be defined as e MeV. The energy 
loss per a nucleon pair is e/{Ap—l) MeV on average. 
Thus, the kinetic-energy loss in the splitting process can 
be written as, 

eAfiA^ - At) , , 

= (6) 

In the case of non-relativistic beam energy, the energy 
conservation between projectile and two pieces gives a 
new formula as follows: 

y_l ^ L £Af{Ap-Af) 

Vp Y ApEpiAp - 1) • ^ ^ 

If we select the energy loss parameter of £=8 MeV, 
this formula corresponds to Borrel's for the case of one- 
nucleon removal. However, the value of e may not be 
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FIG. 7: Peak value, Pq, of fragment momentum distribution 
for HE-component and LE-component. 



8 MeV. In Fig. [6l the Pq values predicted by the for- 
mula were compared with the present data. The solid 
curve (this work) was drawn with e=12 MeV to fit the 
data. The experimental data may support that e=12 
MeV is better than 8 MeV. 



2. Acceleration phenomenon 

The acceleration effect cannot be explained by a frag- 
mentation model based on the abrasion-ablation picture, 
in which the projectile always loses the kinetic energy in 
the laboratory system for abrasion process of nucleons. 

The acceleration phenomenon may be peculiar to the 
^Be^^Ar, X) reaction, namely inverse kinematics. No 
acceleration effect was found in the Ta target data. By 
comparison between both targets, we found two partic- 
ular features of momentum distributions for light frag- 
ments produced in the Ar-fBe reaction. One feature is 
the acceleration phenomenon. The other feature is the 
existence of LE-component. Both phenomena coincide in 
our data. Therefore, the acceleration phenomenon may 
be related to the LE-component. 

F. Auger et al. have suggested that the light frag- 
ments originate in division of a composite system. In 
their experiment of *^Kr+C and ^^Ki+Al collisions, the 
energy of the Kr beam was 35A MeV, and relatively lower 
than our experiment. As a result of their experiment, 
two components of fission-like fragments were found in 
the velocity distribution of fragments with each mass 
number. The two components may be generated from a 
highly excited compound system since the average speed 
of them was equal to the center-of-mass (CM) speed of 
the incident system. In short, it is possible that the 
LE-component and the HE-component come from one 
source. 

As shown in Fig. [71 the velocities of the LE-component 
observed in this work are not below but just on the CM 
velocity. In addition, the mean energy of HE- and LE- 
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component is around 80^ MeV. Thus, the observed phe- 
nomenon is different from the work of F. Auger et al. 

Then, we suggest two source model that the LE- 
component and the accelerated HE-component come 
from two sources produced in different reactions, respec- 
tively. One source is an excited compound system from 
projectile and target nuclei running on the CM velocity. 
The source emits light particles homogeneously in angu- 
lar space and becomes a fragment in the LE-component. 
On the other hand, the other source is a hot projec- 
tile. The hot projectile is formed by a strong impact 
of target nucleus to convert the kinetic energy of pro- 
jectile nucleus to the internal energy for the excitation 
of projectile nucleus. The hot projectile decays to two 
components again, which consists of a faster component 
and a slower component than the projectile velocity, like 
the one-source model. This two-source model predicts a 
third component in the momentum distribution of light 
fragments. Such third component was not observed in 
our data. After all, the origin of the acceleration phe- 
nomenon has not been understood so far. 

Neither the acceleration effect nor the LE-component 
were found in the Ta target data. The different situa- 
tion from the Be target comes from the large impact pa- 
rameter. For a fragment mass, the impact parameter in 
the Ar-|-Ta reaction is larger than in the Ar-fBe reaction 
when the fragment is formed by the geometrical cut. The 
hot compound system is not produced in the peripheral 
collision with the large impact parameter. Therefore, the 
exotic phenomena may not be observed in the Ta target 
data. 



B. High momentum side width 

The high momentum side widths ur as a function of 
fragment mass are shown in Fig. [5Ja) and Fig. (H^a), for 
the Be and Ta targets, respectively. The high momentum 
side widths of all the observed fragments are compared 
to the formulation by Goldhaber as follows. 
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FIG. 8: Momentum widths at high momentum side (Be tar- 
get); (a) for OH and (b) for reduced widths according to the 
Goldhaber model. 



main purpose, we do not follow up the deviation further. 
These results are in good agreement with high-energy 
experiments [l^ . 

At relativistic energies, the reduced width cto is inde- 
pendent of the primary beam energy. At lower energies, 
the reduction of ctq has been observed The reduc- 

tion mechanism has been argued by several theoretical 
works, where the reason is for example due to the Pauli 
blocking. Due to the effect, the ctq has the energy depen- 
dence at 30-40A MeV [1^, [S^l , and becomes constant up 
to 90A MeV. The fact that the measured ctq is the same 
as high-energy one is consistent with this picture. 



where Ap j is the mass number of the projectile and frag- 
ment, respectively. The solid curves are drawn with the 
reduced width (To=90 MeV/c from experimental results 
at relativistic energies [l^l- The deduced values of cto 
for the fragments in the mass range 9 to 36 are shown 
in Fig. [Hl^b) and Fig. [Hl^b). The dashed lines denote the 
mean value of ctq =93.5±2.6stat±7.5ij^s for the Be target, 
and 97.4±1.8stai±7.8sys MeV/c for the Ta target, respec- 
tively. No significant difference between both targets was 
observed in the high momentum side widths. The light 
fragments A/<13 produced by using the Be target show 
slightly a deviation from the Goldhaber model. The de- 
viation is originating from the LE-component and the fit- 
ting function to be used. To avoid a digression from the 



C. Low momentum side width 

The results of momentum width at low momentum 
side are shown in Fig. 1101 The widths of low momen- 
tum side ctl are plotted as a function of fragment mass 
with the statistical errors. We compared the results with 
the high momentum side width an- Instead of showing 
individual data of cjh, the dashed curve calculated with 
the Goldhaber model is presented. The systematic error 
of (Tl is estimated to be 8%, which is not shown in the 
figure. 

In the Be-target data, each isotope chain has a 
mountain-style structure (solid curves). On the other 
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FIG. 9: Momentum widths at high momentum side (Ta tar- 
get); (a) for (TH and (b) for reduced widths according to the 
Goldhaber modeL 



hand, the low momentum widths of the Ta data may 
have no such a structure. One can understand the de- 
viation as the target thickness effect. Another reason 
for the mountain-style structure should be noted. The 
momentum distribution of light fragments has the LE- 
component at the low momentum side only for the ^Be 
target data. In the fitting procedure by use of the asym- 
metric Gaussian-like function, the LE-component may 
also affect the results. Therefore, the mountain-style 
structure of ctl for fragment mass should be ignored as a 
systematic error. 

As seen in Fig.[TOl the measured widths of ctl are twice 
wider than (Jn approximately. Obviously, the large width 
cannot be explained by the Goldhaber model. The 
may have a linear dependence as a function of mass loss 
AA — Ap — Af, which is very different from the parabolic 
feature of an- No target dependence of ctl between the 
Be and Ta targets was found in our data. As any models 
for the low momentum tail have not been established, it 
is difficult to discuss the low momentum tail only from 
the systematics of a^- Yet we should comment that the 
low momentum tail may not depend on the target. 

In addition, another interesting point was found in the 
result. At the limit AA— >0, the ul may not converge on 
0, but 300--400 MeV/c. This feature is different from the 
Goldhaber formula. Even if the fragmentation reaction is 
the dominant process for production of fragments, these 
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FIG. 10: Momentum widths at low momentum side by using 
(a) the Be target and (b) the Ta target. The solid lines in 
both figures are the hnear fitting of the Ta target data. The 
experimental results (a^) and the fitting hues are located far 
from a prediction of the Goldhaber model. The isotope chains 
of experimental data are also drawn as solid lines. 



new observations may lead to an additional reaction pro- 
cess necessary for intermediate-energy reactions. 

In Fig. [TOl a fitting line was first obtained from the Ta 
data. The fitting lines are drawn in both figures. The 
experimental results (ctl) and the fitting lines are located 
far from a prediction of the Goldhaber model. 

What is the origin of the large width of low momen- 
tum side? The large width is indeed produced by an 
energy-loss process in the nuclear reaction. The energy- 
loss process for may be different from the "pure" frag- 
mentation process, because the width of low momentum 
side strongly depends on the beam energy At relativis- 
tic energies, the momentum distribution becomes sym- 
metric with ctl = (7h- For a low energy beam at 30 A 
MeV, the large tail of low momentum side appears obvi- 
ously [IHl. Macroscopic friction process of nuclear dissi- 
pation is not the satisfactory mechanism for because 
the momentum peak Pq, which has no significant energy 
dependence, may also be changed at the same time. 

The energy-loss process to produce the width of low 
momentum side may be explained by nucleon exchange 
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reaction between target and projectile f36|. For instance, 
transfer mechanism adds one nucleon to the projectile 
or fragment, the energy per nucleon E may be changed 
as E+in — E ■ A/{A + 1), where A is the mass number 
before transfer reaction. Since the velocity of projectile 
is reduced by the nucleon exchange reaction, the nucleon 
transfer may contribute the low momentum tail. When 
the projectile gives a nucleon to the target nuclei, the ve- 
locity of the remained nucleons in the projectile does not 
change significantly. Like the effect of momentum peak 
shift, the change of potential energy in projectile affects 
slightly the velocity of projectile nucleus. On the other 
hand, when the projectile picks up a nucleon from the 
target nuclei, the projectile velocity should be deduced 
because nucleons in target nuclei are much slow on av- 
erage at the laboratory system. Thus, we suppose the 
nucleon exchange reaction as an origin of the low mo- 
mentum tail. 

Assumed that the probability of nucleon exchange is 
described as the Poisson distribution and the average of 
number of transfer nucleons in a reaction is quite small, 
the transfer process does not contribute the peak shift 
of fragment-momentum distribution but the large width 
of low momentum tail. If the probability of transfer 
process is small on average, the peak shift of fragment- 
momentum distribution, which is related to the case of 
no nucleon transfer, does not suffer from the transfer pro- 
cess to be independent of the beam energy. On the other 
hand, the low momentum width is sensitive to the trans- 
fer probability. 

We next discuss the linear dependence of ctl as a func- 
tion of Is. A. First, we try to explain this dependence in 
terms of the surface abraded. If the transfer probability 
is proportional to area of the surface, the probability may 
be described as a symmetric function with respect to the 
half of projectile mass. However, the observed behavior 
has a linear dependence. So, the abraded surface may 
not be directly related to the transfer process. 

Secondly, we think about the overlap volume of pro- 
jectile and target nuclei, i.e., total number of nucleons 
in the participant region. The region has almost linear 
dependence as a function oi Af. Thus, the linear depen- 
dence of ctl may be related to the volume of the overlap 
region. 



D. Transfer-like fragmentation 

The transfer-like fragments of 3'^Al,3^'^^Si,^^'^^P, and 
37-40j^ were observed at this experiment. These frag- 
ments have more neutrons or protons than the projectile 
nucleus. The fragments cannot be produced with projec- 
tile fragmentation reactions. Neutron and proton pick-up 
processes are necessary for production of the fragments. 

Figure [TT] shows the momentum distribution of ^^Al 
acquired by using the Be target. This isotope produc- 
tion needs at least one neutron pick-up process from tar- 
get nucleus. We made an attempt to fit the momen- 
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FIG. 11: Momentum distribution for '^^Al. This isotope can- 
not be produced with only the nuclear fragmentation process. 
One neutron should be picked up from target. The measured 
peak momentum corresponds to 83. MeV. The measured 
momentum peak shows a larger shift toward low momentum 
side, compared with that of the same mass number, ~0.9A 
MeV. 



turn distribution with the asymmetric Gaussian function. 
Due to the lack of data at low-momentum tail, we as- 
sumed that the fitting parameters of momentum widths 
were fixed. The an was given from the Goldhaber for- 
mula with (7o=90 MeV/c. The (Jl was assumed to be 
(400 ± QOsys) MeV/c. This value was obtained from the 
systematics as shown in Fig. [TUl 

The fitting result is shown as a dotted curve in Fig. [11] 
The peak of the momentum distribution corresponds to 
83.9^ MeV. As the primary beam energy was measured 
as 87. 5A MeV, the momentum peak shift was obtained 
as (3.6±1.2) A MeV. 

The measured momentum peak shows clearly a larger 
shift toward low momentum side, compared with that of 
the same mass number, '^O.OA MeV. The large momen- 
tum shift of transfer-like fragments was also observed in 
the fragments of s'^'^^Oj^ requiring to proton pick-up pro- 
cess in production. 



E. Target dependence of cross sections 

The target dependence of cross sections was investi- 
gated with the results of our experiments. We have ob- 
tained the cross sections in wide range of fragment charge 
for each fragment mass with small statistical and system- 
atic errors, and for the same projectile (^"Ar) with two 
sets of targets (^Be and ^^-'^Ta). So we can investigate 
the validity of factorization for fragmentation reactions 
at intermediate energies. 

Figure [T^] shows the ratios of cross sections for a frag- 
ment in Ar+Ta reactions to those in Ar-|-Be reactions. 
The cross sections are normalized with the experimen- 
tal mass yield Y{Af ) to eliminate the target size effect. 
The mass yields were obtained from the sum of fragment 
cross sections with the same mass number. The ratio 
is shown as a function of charge difference between the 
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FIG. 12: Ratio of the production cross sections for each 
fragment produced with Be and Ta targets. The sohd hues 
are drawn for the same mass number of nuclei, i.e. isobars. 
The Zp is the /3-stable charge for each isobar. The cross- 
section ratio of (TTa(A,Z)/(Tse(A,Z) is normalized with the 
mass-yield ratio of yTa(A)/Y'se(A) observed. 



most stable charge Zp{Af) and the fragment charge Zf. 
Thus, the target dependence of projectile fragmentation 
cross sections except the target size effect is represented 
for the wide range of isotopes near and far from the sta- 
bility line. If the factorization is valid for production of 
an isotope, the ratio has no Zp — Zf dependence. 

In Fig. [T^l the ratios near the /3-stability line are con- 
stant. It corresponds that the factorization is valid for 
the production of the nuclei. According to the previ- 
ous work [371] . the factorization is valid for the isotopes 
in Zp — Zf<2. On the other hand, the ratios increases 
when Zp — Zf is increased. This deviation shows that 
the factorization hypothesis is clearly broken down for 
neutron-rich nuclei with (Zp — Zf)>2. 



F. EPAX formula 

The EPAX was developed by using data of the spalla- 
tion reactions and heavy-ion induced fragmentation re- 
actions at several A GeV [3^. The cross section of a 
fragment with mass A and proton number Z produced 
by projectile fragmentation from a beam {Ap, Zp) im- 
pinging on a target (At, Zt) can be written as, 

a{A,Z) = Y{A)-W{A,Z) (9) 
W{A,Z) = n-cxp{-R-\Zprob- zf) (10) 

where Y{A) represents the mass yield which is the sum 
of the isobaric cross sections for fragments with mass 
number A^ and W{A,Z) describes the charge distribu- 
tion which means the cross-section distribution of a given 
fragment mass which has a maximum peak at Zprob- 

The charge dispersion W{A, Z) is described as i?, Zprob 
and U parameter. The most probable charge, Zprob, is 



written as, 

Zprob = Zp + + Am, (11) 

where Zp is the /3-stable charge for a fragment of mass 
number, A [39|, A is a proton-excess between the sta- 
bility line and the most probable line of fragmentation 
reaction, and is the so-called "memory effect", i.e., 
the influence of the projectile N/Z ratio on the fragment 
N/Z ratio. The i? is a function of fragment mass num- 
ber, A, which shows a fragment mass dependence of a 
steepness controlling cross sections of isotopes from the 
stability line to the drip lines. The mass dependence has 
been confirmed from many combinations of projectiles 
and targets. The U parameter is given as a constant for 
all neutron-rich isotopes so far. 

Using new data obtained mainly at the GSI/FRS fa- 
cility, the EPAX was recently modified Q slightly to 
tune the mass yield and the U parameter for proton-rich 
fragments in the vicinity of projectiles. The present pa- 
per focus on the neutron-rich fragments in a wide range. 
Thus, this work gives the analysis and discussion using 
the EPAX formula based on the original EPAX that has 
relatively simple functions. 

The EPAX formula is valid for the "limiting fragmen- 
tation" regime, where the fragmentation process is no 
longer energy dependent. The energy dependence of frag- 
mentation cross sections has been investigated by Silber- 
berg and Tsao [13, lU . We can see the similar energy 
dependence in the total reaction cross section. The total 
reaction cross section has been studied both theoretically 
and experimentally for more than 50 years. A detailed 
list of reference is found in Ref. [i^. Compare the cross 
section of ^^C+^'^C at 90A MeV with 900A MeV, the 
difference of cross sections is about 30% [11] . This value 
is nearly equal to the systematic error of measured cross 
sections in this work. Thus, we can assume that the lim- 
iting fragmentation hypothesis is valid in this work. 

The EPAX formula follows the factorization hypoth- 
esis. The mass yield has a target dependence, however 
it is limited to the target nuclear-size effect. The charge 
dispersion VF(A, Z) is independent of the target nucleus. 

In this work, however, the BOF has been found in the 
production of very neutron-rich nuclei. An investigation 
of the charge distribution is necessary to find appropriate 
description for cross sections from our data. In the charge 
distribution of EPAX, the U parameter is a constant of 
1.65 for neutron-rich side, and has no target dependence. 
The value of U parameter is very sensitive to the produc- 
tion cross sections of the isotopes far from the /3-stability 
line. In the following, we reduce the U parameter for 
each fragment mass with both targets from our data. At 
the same time, we investigate the target dependence of 
Zprob since Zprob is also sensitive to the cross sections of 
isotopes far from the /3-stability line. 
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FIG. 14: Deviation of the most probable charge from the 
EFAX parametrization to experimental data for the produc- 
tion targets of Be and Ta. 



A fitting of the production cross sections was per- 
formed using the EPAX function. The function is rep- 
resented as the product of Y{A) and W{A,Z). The 
charge distribution Z) is characterized with the 

most probable charge Zprot, the slope constant U and 
the width parameter R. The R parameter has slightly 
fragment-mass dependence and an effect to the slope of 
charge distribution. We have two kinds of U parameters 
of Up for proton-rich side and J7„ for neutron-rich side. 
In the fitting procedure, the values of the Up and the 
R parameters are fixed, originally given by the EPAX 
formula. The fitting procedure was performed with the 
experimental data of each fragment mass. First, we ob- 
tain the maxima of the charge distribution Zprof)(EXP) 
as a function of fragment mass number from the experi- 
mental data with both targets. The values are compared 
with the EPAX ones, the /3-stability line, and N/Z ratio 
of projectile nuclei. Next, the U parameters are deduced 
from the data when we assume the same parametrization 
of Zprob for both of the targets. 

We have investigated the target dependence of Zprob- 
The measured cross sections for each fragment mass were 
fitted with Y{A) x W{A, Z). The mass yield 1" is a con- 
stant for a fragment mass. Fitting parameters were U 
parameter in neutron-rich side, Zprot, and the mass yield. 
Figure [13] shows the charge distribution of fragment mass 
29, produced in Ar-|-Ta reaction. The Zprob is obtained 
as 13.47±0.01(siai)±0.04(si/s). Similarly, the Zprob value 
for each fragment mass is deduced from our experimental 
data. 

The experimental Zprot is compared with the value 
from the EPAX formula. Figure [14] shows the deviation 
of Zp™;,(EXP) from Zprot(EPAX) for Ar-|-Be and Ar-fTa 
data. Solid lines are drawn to guide the eye. The EPAX 
reproduces the Zprot very well, especially for the mass 
number between 20 to 35. The deviation of Zprot is less 
than 0.2. Due to the lack of proton-rich side data, we 
obtained only three values of Zprot for the Ar-|-Be data 
in Af=18, 22 and 26. To compare both data of Be and Ta 
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FIG. 15: Peak of charge distributions as a function of frag- 
ment mass for the production targets of (a) Be and (b) Ta. 
The solid curves are the most probable charge Zprob of EPAX. 
The dashed curves are the /3-stable charge Z^. The dashed 
lines are the N/Z ratio of *°Ar projectile and compounds of 
"^Ar+^Be and ''°Ar+^'*^Ta, respectively. 



targets shows no significant target dependence of Zprot- 

Figure [15] shows the most probable charge in N/Z 
unit as a function of fragment mass for the Be and 
Ta targets. The most probable charge Zprot from the 
EPAX is close to the /3-stable charge Zp. The ratio 
Np/Zp of the projectile is 1.22. For the composite sys- 
tem of projectile and target, the ratio is represented as 
{Np + Nt)/{Zp + Zt), of which values are 1.23 and 1.43 
for Ar-I-Be and Ar-fTa, respectively. The Zprof<(EXP) is 
well described as the Zprof)(EPAX) as well as Z^. The 
difference of .Zprof)(EPAX) and is mainly the memory 
effect A„i. The memory effect is not clearly seen in our 
data. 
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No significant difference of Zprob(EXP) between tfie Be 
and Ta targets has been observed, however the BOF has 
been found in the very neutron-rich nuclei. What does 
it mean? It should be noted that the BOF had never 
been found in the projectile fragmentation at other ex- 
periments. This fact is consistent with our observation of 
^prob (EXP) for the Be and Ta targets. To understand the 
charge distribution, we show the Zprob of projectile-like 
fragments produced by the reactions at low and interme- 
diate energies. 

At low energies, the most probable charge shows 
the existence of two different reaction mechanisms. 
In Ref. CI isotopes were produced in 7A-MeV 

''°Ar-|-^°Ni reaction and the contour plots of the CI iso- 
tope yield were drawn to the fragment mass and the ki- 
netic energy. We clearly see two components, one cor- 
responding to quasi elastic reactions centered at a high 
energy and a mass of 39. The second one is centered 
around a mass of 36 and a small kinetic energy. The 
most probable charge obtained from the mass number is 
near the N/Z of both projectile and composite system 
which corresponds to the quasi elastic and deep inelastic 
reactions, respectively [IBJ. 

An intermediate composite system has been shown as 
the result of a complete damping of the relative motion 
between the projectile and target nuclei. The projectile- 
like fragments are produced via binary nuclear system 
that collective effects dominate. On the other hand, at 
high energies, the process is dominated by individual nu- 
cleonic collisions that described as participant-spectator 
models. The values of most probable charge distinguish 
the reaction mechanisms to produce the projectile-like 
fragments. 

At intermediate energies, D.Guerreau et al. [i^ re- 
ported the observation of the systematic shift of isotope 
distributions between two targets. The isotope distri- 
butions of fragment yield were measured in *°Ar-fNi 
and ^°Ar-|-Au reactions at HA MeV. The systematic 
shift of the isotope distributions for a given element of 
Si in ^°Ar-|-Au reaction was observed about 0.3 mass 
unit towards the neutron rich side. However, their re- 
sult is different from our data. The fitting results of our 
charge distributions have shown no significant difference 
of Zprof)(EXP) between the Be and Ta targets. 

We have already found that the factorization assump- 
tion is invalid for production of neutron-rich nuclei. The 
flexibility of the charge distribution except the most 
probable charge is the slope parameter. Thus, we would 
seek for the origin of target dependence in the slope con- 
stant U. Due to the BOF, the U may change from a 
simple constant to a complex parameter which depends 
on target nuclei. 

As the target dependence of Zprob has not been found, 
we assume now that the Zprob can be described as the 
EPAX formula for both of the targets. We make an at- 
tempt to fit the data with the function of charge distri- 
bution. Fitting parameters were U in neutron-rich side 
and mass yield Y, and fragment mass dependence of U is 
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FIG. 16: U parameters as a function of fragment mass for 
the production targets of Be and Ta. On the assumption that 
the value of U is same for. every group of four sets of Af 
data, the fitting results are shown with the systematic error. 
The dashed lines are the values of (7=1.65 from the EPAX. 



tried to be obtained. The fitting procedure is performed 
for the data of each fragment-mass group, where every 4 
mass units are combined to avoid poor statistics. 

Figure [12] shows the U parameters as a function of 
fragment mass for both targets. We plot the values with 
the systematic errors. The slope parameter U of EPAX 
is a constant of 1.65 (dashed lines). The average of U is 
1.62 for the Ta target. The smaller U gives larger values 
of production cross sections for neutron-rich nuclei. On 
the other hand, the U parameter of Be target is larger 
than 1.65 and shows fragment mass dependence. Our 
data demonstrate the BOF for very neutron-rich nuclei 
as the target dependence of U parameter. 



2. Predictive power of new parametrization 

We have obtained the modified EPAX formula for 
the nuclear fragmentation at an intermediate energy for 
both Be and Ta targets. To confirm the validity of 
our parametrization, we show an example of the pro- 
duction cross section predicted for extremely neutron- 
rich nuclei. Figure fTTl shows the predictive power of the 
new parameterization. The dashed and solid curves are 
the charge distribution of mass-24 isobars produced in 
""^Ar-l-^Be and '^''Ar-f ^^^Ta reactions, respectively. The 
solid box is the cross section of fragmentation channel for 
^*^Ta(^Ar,^^0). This value acquired in another experi- 
ment [47| is in agreement with our new parametrization 
clearly. 
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FIG. 17: Predictive power of cross sections of the new 
parametrization. The solid box is the data of ^'*0 with a 
tantalum target. 



G. Mechanism of Prefragment Production 

We have investigated the nuclear fragmentation reac- 
tion in nucleus-nucleus collisions at intermediate ener- 
gies. The momentum distributions of projectile-like frag- 
ments have been measured for various isotopes includ- 
ing very neutron-rich nuclei. The charge distributions 
of fragment cross sections acquired with the Be and Ta 
targets have revealed the BOF for very neutron-rich nu- 
clei. In order to explain the origin of BOF, we discuss the 
mechanism of prefragment production in terms of the AA 
model. In the AA model, the production cross section of 
a final fragment is determined by the charge distribution 
and the excitation energy of prefragments. 

First, we discuss the charge distribution of prefrag- 
ments on the basis of AA model. The abrasion model is 
a macroscopic description which gives a picture of a clean 
cut of the projectile nucleus by the target nucleus [is! ]. 
The concept of geometrical separation in an overlap or 
participant zone does not specify its proton-to-ncutron 
ratio. A few different methods to calculate the charge 
distribution of the prefragments are so far pro pos ed, a 
correlation model (49| , hypergeometrical model [50| , and 
GDR model [Slj . After all, the charge distribution of pre- 
fragments from these methods depends on the projectile 
nuclei but has no target dependence. 

Next, we discuss the excitation energy of prefragments. 
In an early model [s^l, the average excitation energy of 
a prefragment after abrasion was estimated from the ad- 
ditional surface energy due to the excess surface area. 
Several studies have been made on the source of the ex- 
citation energy [1, [sl, [13]. According to the researches, 
the surface energy is not the main source of the excitation 
energy, but, one of the energy sources for prefragment 
excitation. At least, the target dependence of the excita- 
tion energy can be seen in the surface energy. Figure [T51 
shows a schematic view of prefragment production. The 
prefragments of the same mass number A f ' are produced 



FIG. 18: A schematic drawing of prefragment production. 
The prefragments of the same mass number {Af') are pro- 
duced by projectile fragmentation reaction with target nu- 
clei of dilferent nuclear radius. For the large target nuclei, 
the large impact parameter leads to (A) gentle reaction, and 
(B) violent reaction is expected for the small impact param- 
eter with the small target nuclei. 



by the projectile fragmentation reaction with target nu- 
clei of different nuclear radius. If a target nucleus has a 
large size, the impact parameter is also large to produce 
the prefragment (case A). On the other hand, the small 
size of a target nucleus caused the violent collision to 
form the prefragment with the same mass number (case 
B). In the geometrical model, we calculated the surface 
energy of prefragments in Ar+Be and Ar-|-Ta reactions. 
As the result, large difference of Eg is presented for both 
of the targets. As for the prefragments of AA = 20, the 
values of Eg are approximately 90 MeV and 10 MeV for 
the Be and Ta targets, respectively. The target depen- 
dence of prefragment excitation energies stems from the 
additional surface energies in the AA model. 

We observed no significant target dependence of mo- 
mentum peak shift in the data for the projectile- like frag- 
ments of Af >20. Assumed that nucleons are emitted 
homogeneously from the hot prefragment in the ablation 
process, the velocity of a prefragment does not change on 
average when it becomes the final fragment. The results 
of momentum peak shift show that the energy losses in 
the fragmentation process for the fragment of = 20 
were 60~80 MeV for both of the targets. The energy 
losses are not correlated to the additional surface energies 
on the calculation. Thus, we cannot deny the possibil- 
ity that the BOF originates from the target dependence 
of the charge distributions for prefragments except the 
excitation energies. 

To apply this idea to the BOF problem, we suggest two 
mechanisms of a fluctuation giving for the charge distri- 
butions of prefragments. One is the nuclcon exchange 
process during the fragmentation process, and another 
is the effect from the difference of nucleon-nucleon total 
cross sections {crpp,<Tnn and (Tp„). 

We observed the nuclcon exchange process during the 
fragmentation process in our data. The fitting results of 
Pq and an carry the information at the abrasion stage 
and supporting that "pure" projectile fragmentation is 
a dominant process at intermediate energies. As for cji, 
the larger values oi <tl than cth may be related to the 
energy loss due to the transfer process. Unlike the be- 
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havior of (Th , ctl has a linear dependence as a function of 
fragment mass. The linear dependence may be connected 
from both projectile and target nuclei. So, the source of 
transferred nucleons can be assumed to be the participant 
region. Since the Ta target nucleus is more neutron-rich 
than the Be, the probability of neutron transfer should be 
large. The nucleon exchange reaction between projectile 
and target nuclei during the abrasion process may break 
the factorization of fragment production cross sections. 

The difference of the total cross sections between 
proton-to-proton {(Jpp) and proton-to-neutron (crp„) may 
bring the fluctuation of prefragment charge distributions. 
As a fact known well, the cr„p is larger than the (t„„ 
and CTpp when the nucleon energy is less than 500 MeV. 
In projectile fragmentation process, a neutron-rich tar- 
get nucleus easily knocks out protons in the projectile 
nucleus, so that the production of neutron-rich prefrag- 
ments is promoted. 

The discussions so far are qualitative. However, they 
gives several predictions. First, if the nucleon exchange 
process is just the reason of BOF, the probability of nu- 
cleon exchange reaction becomes large at low energies. 
The nuclear fragmentation experiments with lower inci- 
dent beam energies than at our experiment may reveal 
the large BOF effect for very neutron-rich nuclei. Sec- 
ondly, if the difference of N-N reaction cross sections 
causes the BOF, the nuclear fragmentation experiment 
with 100-800A McV beams should be carried out to con- 
firm the BOF for very neutron-rich nuclei. Since the cross 
section curves of cFpp and fTp„ cross at 500^ MeV, the 
difference of U-paramctcr between the Be and Ta targets 
might change to be reversed. The further investigations 
of fragmentation cross sections at around 500A MeV are 
of great interest. 

V. SUMMARY AND CONCLUSION 

The projectile fragmentation reactions at intermediate 
energies have been investigated using a 90-94^4 MeV 



^"Ar beam at RIKEN-RIPS. We paid our attention to 
the target dependence. Measurement of longitudinal mo- 
mentum distributions of projectile-like fragments within 
a wide range for fragment mass and its charge including 
very neutron-rich nuclei has been performed with ^Be 
and ^®^Ta targets. From the momentum distribution of 
fragments, a parabolic mass dependence of momentum 
peak shift was observed in the data of both targets, and 
a phenomenon of fragment acceleration was observed 
only in the Be-target data. A linear dependence of the 
low momentum tail as a function of removed nucleons 
was found for both targets. As a possible origin of the 
low momentum tail, we suppose the nucleon exchange 
reaction. We observed large target dependence of the 
cross sections to produce very neutron-rich nuclei with 
{Zp — Zf)>2. The deviation shows that the production 
cross sections of very neutron-rich nuclei far from the 
stability line do not factorize. The reaction mechanism 
was discussed, however the discussions so far are quali- 
tative. It is much interesting to investigate the target 
dependence of the cross sections at higher energies. 
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